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Pseudogap effects induced by resonant pair scattering 
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We demonstrate how resonant pair scattering of correlated electrons above Tc can give rise to pseu- 
dogap behavior. This resonance in the scattering T-matrix appears for superconducting interactions 
of intermediate strength, within the framework of a simple fermionic model. It is associated with 
a splitting of the single peak in the spectral function into a pair of peaks separated by an energy 
gap. Our physical picture is contrasted with that derived from other T-matrix schemes, with super- 
conducting fluctuation effects, and with preformed pair (boson-fermion) models. Implications for 
photoemission and tunneling experiments in the cuprates are discussed. 



PACS numbers: 74.20.Mn, 74.25.-q, 74.25. Fy, 74.25.Nf, 74.72. 



cond-mat/9705144 



It has become clear in recent years that the pres- 
ence of a pseudogap above the superconducting transition 
temperature, Tc, is a robust feature of the underdoped 
cuprates. This phenomenon is manifested in thermody- 
namic |l| , magnetic [|j , and angle-resolved photoemission 
spectroscopy (ARPES) data g. These ARPES experi- 
ments, which have established the presence of a Luttinger 
volume Fermi surface, place important constraints on any 
pseudogap scenario: they indicate that the pseudogap 
appears directly in the spectral function and its magni- 
tude and symmetry seem to evolve smoothly into that 
of the superconducting state. Furthermore, the minimum 
gap points in the pseudogap regime retrace the normal 
state Fermi surface Q. 

A variety of theoretical scenarios have been proposed, 
for the origin of the pseudogap. Quantum Monte Carlo 
simulation studies have been carried out on both positive 
and negative U Hubbard models [||. Alternative models 
relate the pseudogap to either magnetic pairing of spins 
RVB-like pairing of chargeless spinous 0, or pre- 
cursor superconductivity effects The present paper 
addresses this last scenario, in part because of constraints 
from ARPES data and in part, because the cuprates are 
short coherence length, quasi-two dimensional supercon- 
ductors, with anomalously low plasma frequencies 
They are, therefore, expected to exhibit important devi- 
ations from an abrupt, BCS-like transition. 

In our physical picture, we associate an important com- 
ponent of the cuprate pseudogap with resonant scatter- 
ing between electrons of opposite spin and small total 
momentum. This resonance arises in the presence of 
intermediate coupling and a sizeable Fermi surface. A 
depression in the density of states occurs because states 
near this Fermi surface are unavailable for electrons in the 
Fermi sea to scatter into; such states are otherwise occu- 
pied by relatively long lived (metastable) electron pairs. 
The related suppression in the spectral weight differs 
from that derived from conventional low frequency and 
long wavelength fluctuation effects In the present 

case it is the strength of the attractive interaction, rather 
than the critical slowing down (in proximity to Tc), which 



leads to the long-lived pair states. It should be noted that 
our resonant scattering approach is to be distinguished 
from previous precursor superconductivity models asso- 
ciated with either preformed pairs or dynamic phase 
fluctuations 

Our starting point is a scheme which connects the 
strong coupling, short coherence length description of 
superconductivity formulated by Leggett and Nozieres 
and Schmitt-Rink [|l2| with a well established T-matrix 
formalism designed to treat normal state fluctuation ef- 
fects in conventional superconductors We consider 
a generic model Hamiltonian 
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where cj^^ creates a particle in the momentum state k 
with spin a, and = kF' /2m — n (we take h = ks 
Here I/k,k' = g<^k¥'k', where (^k = (1 + fc^/^o) 
and g < is the coupling strength. While we consider 
the s-wave symmetry case, d-wave symmetry can be 
readily introduced via (pk {coskx — cos ky). In the 
scheme of Nozieres and Schmitt-Rink the transition tem- 
perature, Tc, must be obtained in combination with the 
chemical potential, fi, by use of the Thouless criterion, 
^q^o(^ = 0) = (see below), and the usual equation for 
particle number. These authors have demonstrated that, 
when the parameter g is varied, the appropriate coupled 
equations lead to an interpolation scheme which contains 
the BCS limit for small g/gc , where /i ~ Ep, and that 
of Bose-Einstein condensation in strong coupling, where 
fi becomes negative. Here gc — —An/ mho is the criti- 
cal value of the coupling at which a bound state of the 
isolated pair first appears. 

In this paper, we extend this formulation so as to pro- 
vide a basis for computing the spectral function and den- 
sity of states and to simultaneously incorporate appropri- 
ate conservation laws |Q. To this end, we calculate the 
single particle self-energy, Ek((^), and the T-matrix as 
shown diagrammatically in Fig. ^ [|l5|. Analytically, the 
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FIG. 1. Diagrams for coupled E, T, in full scheme (a) and 
lowest order conserving scheme (b) used here. Fig. 1(c) rep- 
resents a rewriting of T in Fig. 1(a). 



self-energy corresponds to 



and the T-matrix, which may be written in the Dyson 
form shown in Fig. |l](c), is given by 



Tq(ir2„) 



p,Ci 



(3) 



Here G^\iQ) = ^''^"'('O) - Sk«/) and G^l\i<:i) = 
{iC,i — €k)^^ ■ Finally, the spectral function is defined as 
^k(w) — — TT^^Im [G'k(«0 ^ 1^ + iO)]; this leads to the 
density of states, N{ijj) — X]k^k(w). 

The choice of diagrams to include in a T-matrix scheme 
has been extensively discussed in the literature [p"5|-pT[. 
Our asymmetric choice [Fig. |l|(a)] - in which the T- 
matrix contains one self-energy renormalized and one 
"bare" propagator - builds on the early work of Kadanoff 
and Martin |16| ]. When this diagrammatic scheme was 
applied to conventional superconducting fluctuation ef- 
fects it was shown ||l5[ that the fluctuation gap above 
Tc smoothly evolved into the superconducting gap be- 
low Tc- Moreover, this approach is known []l5| , p^ to re- 
produce the conventional BCS theory in the appropriate 
weak coupling limit. Direct connection can be made to 
the related theories of Marcelja jl^ and Yamada and col- 
laborators jl^ if the full line in the "box" (the pair fluc- 
tuation self-energy) of Fig. ^(c) is replaced by a noninter- 
acting line. Finally, it is straightforward to demonstrate 
using the more general criteria introduced by Kadanoff 
and Baym that this theory preserves all conservation 
laws. 

It has been shown that the results of the full "mode- 
coupling" scheme of Fig. |l|a are qualitatively captured by 
the lowest order conserving approximation of Fig. 0(b) 
[ p^ ; this approximation, nevertheless, goes beyond the 
original work of Nozieres and Schmitt-Rink. Unlike other 
T-matrix approaches, where higher order self-consistency 
effects tend to diminish leading order features |Q, it is 




1.0 



0.0 



1.0 




-1.0 



0.0 



1.0 



FIG. 2. A]^{uj) vs Lu for weak {g/gc ~ 0.6) in (a) and inter- 
mediate {g/gc ~ 1.0) coupling in (b). T/Tc varies from 1.0 to 
1.1. Insets plot Re [rqllQ(r2)] (solid lines, for same T/Tc, as 
in main figure), and Im [Tqll(,(n)] (dashed lines, at T = Tc). 



found p8[ | that within the present framework the inclu- 
sion of "mode-coupling" effects amplify these first order 
(pseudogap) features. For simplicity, we, therefore, fo- 
cus on the lowest order approximation. The discussion 
of feedback effects is deferred to a future publication. 

In our scenario the physical process which generates 
the pseudogap is resonant pair scattering (above Tc), aris- 
ing from the condition that the real part of the inverse 
T-matrix, Re[Tq"^(r2 = flq)] = 0, when the imaginary 
part, Im [r^^(riq)], is sufficiently small. This resonant 
behavior is manifested as a sharp peak in Im[rq(rj)]. 
This peak is in turn reflected in the electronic self-energy 
and the spectral function. The pair resonance condi- 
tion is illustrated in the insets of Fig. ^, where the be- 
havior of Re [rql^Q(ri)], as a function of frequency is 
contrasted for weak {g/gc < 1) (2a) and intermediate 
{g/9c ~ 1) (2b) couplings. Each series of curves corre- 
sponds to varying temperature. The dashed lines indi- 
cate the form of Im[r^^(ri)] at Tc in each of the two 
cases. The critical value gc establishes the approximate 
dividing point between resonant and non-resonant scat- 
tering. As can be seen, there is a finite frequency zero 
crossing of Re [r^^(r2)] for T > Tc, corresponding to res- 
onant scattering, only in the stronger coupling limit. The 
resonance energy increases as a function of temperature 
T and q until it disappears at a cross-over wave vector 
q* or temperature T*. 
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The associated spectral functions Ak(u;) for each of the 
two cases considered in the insets are numerically com- 
puted from the self-energy using Eq. (||) and plotted for 
the case k = kp in the main portion of Fig. ^ as a func- 
tion of Lu, for varying T. (Throughout the unit of energy 
is Ep.) Although the numerical integrations involved are 
computationally intensive, the integrated spectral weight 
is unity to several significant digits for each spectral curve 
presented. In the stronger coupling limit and at suffi- 
ciently low T (Fig. ||(b)), the two peaked structure char- 
acteristic of a pseudogap appears and becomes more pro- 
nounced with larger g/gc- In the more weakly coupled 
limit {g/gc = 0.6), the single peak behavior character- 
istic of a normal Fermi liquid is recovered, as shown in 
Fig. 11(a). In general, the two peaked structure correlates 
with the presence of a resonance in the T-matrix. For g 
slightly greater than gc, the two maxima are resolvable 
up to T* of the order of several Tc. 

An intuitive understanding of the splitting of the spec- 
tral peak into a pair of asymmetrically broadened peaks 
may be gained by examining the imaginary part of the 
self-energy. On the real frequency axis {iQ uj + iS), 
Im [I]k('^)] is given by 

Im[l]k(w)] = -^^k-q/2lm[Tq(w-(-eq_k)] X 
q 

[/(fq-k) + n{u; + eq_k)] • (4) 

where f{x),n{x) — {e^^ ± 1)^^ For intermediate 

coupling strengths, a resonance condition leads to a peak 
in Im [Tq(r2)] at small frequencies and momenta, which 
in turn yields a maximum in — Im [Ek('-^)] at u; -I- ek ~ 
(see the inset of Fig. ||(b)). The frequency weight under 
this peak is written as 7r|Ap(y9^, where |A| can be viewed 
as the pseudogap energy. This peak in — Im[I]k(w)] im- 
plies - via the Kramers-Kronig relation - a corresponding 
resonance structure in Re [I]k('^)] at the same frequency 
u; « — Ek- In this way Ak(i-Li) acquires two peaks sepa- 
rated by 2|A|(pk with 

/-\-oo 1 (-) 
^nmm[T^m]. (5) 
^ -°° 

The asymmetric broadening of the two spectral peaks 
is a generic feature of our results and is dvie to the in- 
teraction of correlated pairs with the Fermi sea. This 
asymmetry, which is contained in Eq. ^, reflects that in 
Im[Tq(ri)], as a function of fl. In Fig. || we plot the 
momentum dependence of the spectral function slightly 
above Tc for weak (3a) and intermediate (3b) coupling, 
along with typical self-energies shown in the insets. The 
former case shows the single peak structure which evolves 
with k in a fashion characteristic of a finite temperature 
Fermi liquid p^ . In the stronger coupling limit (Fig. ^) 
the spectral weight shifts from the negative to the pos- 
itive frequency peak as the momentum vector k passes 
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FIG. 3. A]^{uj) vs Lu for (a) weak [g — O.Sgc) and (b) inter- 
mediate {g — 1.2gc) coupling, from k < hp (bottom curve) 
to k > kp (top curve). The spectra in (b) show signatures of 
normal state "particle-hole mixing" . In the insets are plotted 
the corresponding Im [Ek(ii')] and Re [Ek(i^)] for k = kp- 

through the Fermi surface. Close to the Fermi momen- 
tum the peaks disperse roughly as = ±\/ e^. + |Ap(^j^. 
This dispersion provides a predictive signature for future 
ARPES measurements, within the precursor supercon- 
ductivity scenario. Indeed, this behavior is reminiscent 
of the particle-hole mixing found in photoemission mea- 
surements on the superconducting state [ p3|p4| ] . 

Finally, the density of states, N{uj), is plotted in Fig. ^ 
as a function of energy. This quantity may be directly 
related to tunneling as well as to thermodynamic mea- 
surements in the pseudogap regime. The asymmetry in 
the curves reflects, in part, the asymmetry of the spectral 
functions seen in Fig. || and ^. For clarity the results are 
represented by subtracting the "normal" state curve, ob- 
tained, for definiteness, in the very weak coupling limit. 
Fig. ^(a) indicates the coupling constant dependence of 
N{uj) and Fig. ^(b) the corresponding temperature de- 
pendence for fixed g. A depression in N{uj) - which in- 
creases with g - develops at smaller couplings, and per- 
sists to higher temperatures, than do pseudogap effects 
in the spectral function (see Fig. ||). 

In summary, we have demonstrated how resonant pair 
scattering above Tc gives rise to a splitting of the spectral 
function. Ay {uj) , as well as a density of states depression. 
Experimental observation of the former is the more sig- 
nificant manifestation of pseudogap behavior, providing 
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FIG. 4. N{ij) vs uj for g/g^ = 
and for T/Tc = 1.0 to 1.2, at g 
density of states is subtracted off in both cases (see text) 



0.8 to 1.1. and T = Tc (a), 
= gc, (b). Tiie normal state 



strong constraints on theoretical models. Our precursor 
superconductivity scenario has predictive signatures: an 
asymmetry in the widths of the two spectral peaks and 
a k-dependent dispersion of the T > Tc spectral func- 
tion, qualitatively similar to that of the BCS state. A 
dx^-y^ symmetry of the normal state gap will arise natu- 
rally in the present scenario, for a d-wave superconduct- 
ing instability. This would be accompanied by a spectral 
peak broadening proportional to (cos — cos ky f . The 
present picture should be differentiated from preformed 
pair models: the correlated pairs of our picture have 
significant spatial extent and fail to obey Bose statis- 
tics. Furthermore, in contrast to the stripe picture of 
Emery and Kivelson, the amplitude and phase of this 
paired state is never established beyond the dimensions 
and lifetime of a single pair. Quasi two dimensionality 
will enhance our pseudogap effects, which should, then, 
become more pronounced as the insulator is approached. 
Magnetic correlations may, also, ultimately play a role 
in the extreme underdoped regime. Nevertheless, short 
coherence lengths and quasi 2d features suggest that pre- 
cursor superconductivity is present to some degree and 
must necessarily be calibrated in order to obtain a full 
understanding of the cuprate pseudogap. 
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